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COVERAGE  STATISTICS  OF  DISTRIBUTED  SENSOR  FIELDS 
WITH  HETEROGENEOUS  RANGE  SENSITIVITY 


1.  INTRODUCTION 


This  report  examines  the  detection  performance  of  a  distributed  sensor  field  composed 
of  a  moderate  to  large  number  of  sensors  having  heterogeneous  range  sensitivities.  The  need 
to  rapidly  deploy  a  large  number  of  sensors  on  short  notice  could  place  limitations  on  sensor 
quality.  In  particular,  the  requirement  that  all  sensors  have  the  same  range  sensitivity  could 
be  economically  prohibitive  for  even  moderately  large  systems.  This  report  examines  the 
coverage  statistics  for  a  system  of  randomly  distributed  sensors  having  heterogeneous  range 
sensitivities.  Throughout,  comments  are  made  on  the  implications  of  the  analytical  results  on 
design  guidance  and  sensor  deployment. 

This  research  differs  from  previous  work  in  this  area  in  that  it  allows  for  the  range  sen¬ 
sitivity  of  the  sensors  to  be  independent;  it  requires  only  the  distributional  characteristics  of 
the  range  sensitivities  to  be  known.*  Other  work  in  this  area  considers  the  coverage  provided 
by  systems  of  sparsely  distributed  sensors,  all  of  which  have  identical  range  sensitivity;  see 
for  example,  Wettergren  [  1  ]  or  Cox  [2].  In  a  recent  article.  Wan  and  Yi  [3]  examine  the  cov¬ 
erage  provided  by  randomly  distributed  wireless  sensor  networks.  However,  their  research 
considers  only  homogeneous  sensor  systems,  that  is,  systems  in  which  all  sensors  have  iden¬ 
tical  range  sensitivity.  This  analysis  differs  also  from  previous  research  in  that  it  considers 
the  problem  in  both  two-  and  three-dimensional  space;  previous  work  in  the  area  considered 
only  the  two-dimensional  problem. 

Section  2  of  this  report  begins  with  a  brief  review  of  temporal  Poisson  point  processes, 
followed  by  a  discussion  of  spatial  Poisson  point  processes.  The  main  result  on  the  cover¬ 
age  of  a  randomly  distributed  sensor  Held  is  given  in  section  3.  Section  4  contains  examples 
demonstrating  how  the  theory  can  be  used  to  measure  sensor  Held  coverage.  Section  5  com¬ 
pares  the  coverage  of  a  sensor  field  with  heterogeneous  range  sensitivities  to  the  coverage  of 
a  sensor  field  in  which  all  sensors  have  identical  sensitivity  range.  Section  6  is  a  simulation 
study  that  verifies  the  theoretical  results.  Finally,  appendix  A  contains  a  conditional  uniformity 
of  sensor  locations  and  appendix  B  contains  a  proof  of  the  main  result. 

*It  is  supposed  that  the  distributional  characteristics  of  the  sensors  will  be  provided  by  the  sensor  manufacturer  or  can  be  determined 
through  sampling. 
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2  SPATIAL  POISSON  PROCESSES 


This  section  commences  with  a  review  of  the  more  familiar  temporal  Poisson  point  pro¬ 
cess,  followed  by  a  description  of  its  spatial  equivalent.  Recall  that  the  usual  (i.e.,  temporal) 
Poisson  point  process  Nt,  with  rate  A,  is  a  stochastic  process  in  which  Nt  is  the  number  of 
occurrences  of  an  event  during  t  units  of  time.  Moreover,  Nt  follows  a  Poisson  distribution 
with  mean  A t.  Finally,  the  event  occurrences  in  disjoint  time  intervals  are  independent.  More 
formally,  the  counting  process  N  =  {Nt,t  >  0}  is  a  Poisson  process  with  rate  A  >  0,  if  it 
possesses  the  following  properties:* 

1 .  N0  =  0, 

2.  P{Nh  =  1}  =  A h  +  o{h), 

3.  P{Nh  >  2}  =  o{h), 

4.  Nt  satisfies  the  stationary  and  independent  increment  property. 

Property  2  says  that,  over  a  very  short  time  period,  the  probability  of  an  event  occurrence 
is  approximately  linear  with  time.  Property  3  says  that  the  probability  of  two  or  more  oc¬ 
currences  over  a  very  short  time  period  is  essentially  zero.  And  according  to  property  4,  the 
process  Nt  is  independent  of  the  process  Nt+h  —  Nt,  i.e.,  the  number  of  occurrences  up  to 
time  t  is  independent  of  the  number  of  occurrences  after  time  t. 

For  a  spatial  Poisson  process,  instead  of  counting  the  number  of  occurrences  over  a 
time  interval,  one  is  interested  in  the  number  of  points  associated  with  an  event  over  a  region 
of  space.  Let  S  be  an  n-dimensional  set  and  suppose  A  C  S,  i.e.,  A  is  a  subset  of  S.  Consider 
points  scattered  randomly  throughout  S  and  let  N(A)  denote  the  number  of  points  from  the 
scattered  set  that  are  contained  in  A.  The  stochastic  process  N(A)  is  called  a  point  process  in 
S.  Depending  on  the  dimension  of  <S,  let  ||A||  denote  the  length,  area,  volume,  etc.,  of  A.  The 
stochastic  process  N  =  {N(A),  A  C  5}  is  an  n-dimensional  Poisson  counting  process  with 
parameter  A  >  0  if 

(a)  N(A)  follows  a  Poisson  distribution  with  mean  A||A||,  and 

(b)  the  number  of  points  occurring  in  disjoint  subsets  of  S  are  mutually  independent. 

Thus,  for  the  spatial  Poisson  process,  the  statements  equivalent  to  properties  1  —  4  above  are 
(here  0  denotes  the  empty  set): 

'The  function  g(h)  is  said  to  be  little-oli  h,  written  g(h)  =  o(h)  if  lim/,  _0  =  0.  Thus,  according  to  properly  3,  as  h  goes  to  zero, 

the  probability  of  two  or  more  occurrences  goes  to  zero  faster  than  any  linear  function. 
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1.  N(0)  =  0, 

2.  P{jV(A)  =  l}  =  A||A||+o(||A||^ 

3.  P{N(A)  >  2}  =  o(||A||), 

4.  N(A)  satisfies  the  stationary  and  independent  increment  property. 

In  particular,  property  (2)  says  that  if  ||A||  (the  volume  of  the  set  A)  is  small,  then  the  proba¬ 
bility  of  one  occurrence  is  approximately  linear  with  respect  to  volume;  property  (3)  says  that 
for  small  volumes,  the  probability  of  two  or  more  occurrences  is  approximately  zero.  Finally, 
property  (4)  says  that  if  A,  C  C  S,  then  N(A)  is  independent  of  N(C  fl  Ac). 

The  definitions  are  summarized  as  follows:  IfN  =  {N(A),A  C  <S}  isan  ^-dimensional 
Poisson  counting  process  with  parameter  A  >  0,  then  the  probability  of  k  occurrences  in  the 
subset  A  is  given  by* 


P{N(A)  =  A:} 


(AH/iiDV^ii 

k\ 


See,  for  example,  Resnick  [4]  or  Kingman  [5]  for  further  information  on  spatial  Poisson 
processes. 


"Recall  that  if  X  is  a  Poisson  random  variable  with  mean  A,  then  P{X  =  k}  =  Xke  x/k\.  According  to  property  (a),  therefore,  one 
can  substitute  N(  A)  for  X  and  A  ||  A  ||  for  A  to  get  the  equivalent  probabilities  for  the  spatial  random  variable. 
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3.  SENSOR  COVERAGE 


This  section  contains  the  main  result  of  this  report,  with  the  remainder  of  the  report 
addressing  the  consequences  of  this  result.  Most  of  the  discussion  concerns  coverage  in  three- 
dimensional  space;  however,  all  arguments  apply  to  two  dimensions  also.  A  two-dimensional 
equivalent  of  the  main  result  is  also  stated.  The  proof  is  technical  and  is  postponed  until 
appendix  B.  A  restatement  of  the  problem  can  be  found  in  chapter  16  of  Karlin  and  Taylor  [6], 
where  instead  of  sensors  and  range  sensitivities,  spherical  centers  and  corresponding  radii  are 
considered. 

Theorem  1:  Consider  a  set  of  omnidirectional  sensors  in  three-dimensional  space  scattered 
throughout  some  region  S  C  R 3.  Suppose  that  the  sensors  are  spatially  distributed  according 
to  a  Poisson  point  process  with  parameter  A.  Suppose  also  that  the  range  sensitivity  of  all 
sensors  is  distributed  according  to  the  cumulative  distribution  function  F(r)  with  density 
F'(r)  =  f(r)  having  finite  third  moment.  Finally,  assume  that  range  sensitivity  is  independent 
of  the  sensor  location.  Then,  the  number  of  sensors  that  will  detect  a  target  located  at  some 
point  x  £  S  is  a  Poisson  random  variable  with  mean 


Remark  1:  Given  an  arbitrary  point  x  6  <S  C  Rn , 


Afce  ^ 

Pjtarget  at  position  x  detected  by  A:  sensors}  =  — — — . 

The  probability  that  the  point  is  not  detected  is  exp  (—A).  Hence,  the  probability  that  at  least 
one  sensor  covers  the  point  is  p  =  1  —  exp  (—A). 

Definition  1:  Throughout  the  remainder  of  this  report,  the  statistic  p  =  1  —  exp  (—A)  will 
be  referred  to  as  the  predicted  coverage  for  the  sensor  network.  The  coverage  statistic  p  is  a 
measure  of  the  probability  that  an  arbitrary  point  is  detected  by  at  least  one  sensor. 

Remark  2:  If  an  arbitrary  point  in  the  volume  is  to  be  detectable  by  at  least  some  fixed  number 
of  sensors,  then  the  result  shows  us  what  intensity  level  (i.e..  A)  is  required  to  achieve  this. 
Since  the  mean  number  of  sensors  that  detect  a  target  at  some  point  in  S  is  |  A7r  r3f  ( r)dr , 

one  can  adjust  A  (by  increasing  the  number  of  sensors  scattered  throughout  S)  to  achieve  the 
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desired  coverage.  The  only  quantity  that  is  required  is  some  estimate  of  /(?■),  the  range 
sensitivity  density. 

Remark  3  [The  case  of  fixed  range  sensitivity]:  The  reference  to  range  sensitivity  distribu¬ 
tion  alludes  to  the  fact  that  the  sensors  are  not  perfect.  Thus,  instead  of  a  few  very  expensive 
sensors  with  fixed  range  sensitivity  rG,  the  network  may  contain  several  less  expensive  sen¬ 
sors  with  variable  range  sensitivities.  Nevertheless,  through  sampling,  one  can  get  an  estimate 
of  the  sensitivity  distribution.  In  the  idealized  case  in  which  all  sensors  have  identical  fixed 
range  sensitivity  r0,  one  obtains 


3.1  NEAREST-NEIGHBOR  PROBLEMS 

Suppose  that  sensors  are  distributed  throughout  some  three-dimensional  space  with  un¬ 
known  or  slowly  changing*  Poisson  parameter  A.  If  the  sensors  are  designed  to  communicate 
with  each  other  so  that  it  is  possible  to  estimate  the  average  distance  between  them  at  any 
time,  then,  as  will  be  shown  presently,  this  result  can  be  used  to  determine  the  current  value  of 
A. 

Proposition:  Consider  omnidirectional  sensors  in  FtA  distributed  according  to  a  Poisson  pro¬ 
cess  with  intensity  parameter  A.  The  distribution  function  Fo(r)  of  the  distance  between  a 
sensor  and  its  nearest  neighbor,  and  the  mean  distance  E(D )  between  nearest  neighbors  are, 
respectively,1' 


FD(r)  =  1  -  exp 


and 


E(D)  = 


r(|) 

\/36A7r 


Proof:  Let  A  be  the  region  of  space  within  a  distance  r  of  the  reference  sensor.  Then, 


Fo(r)  = 


=  1  —  P  {nearest  neighbor  not  within  distance  ?  } 

=  1  —  P  {no  other  sensor  within  distance  r  of  this  one} 

=  1  -exp{-A||A||}  iiMliA! 

f  4  ,} 

=  1  —  exp  <  — A-7rr  > . 


*  For  example,  the  parameter  A  may  change  due  to  ocean  currents, 
t  Recall  that  r(^)  =  tz~1e~tdt. 
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The  nearest-neighbor  density  /D(r)  =  F'D(r)  =  4A7rr2exp  {— A|7t?’3}  .  Therefore, 


3T(|) 


^36Att 

r(|) 


(2) 


v^36A7t 


where  the  identity  T(r  +  1)  =  rT(r)  for  r  >  0  is  used  in  the  last  line. 

Knowing  A  enables  one  to  determine  the  probability  of  detecting  an  object  that  enters 
the  sensor  network  space:  From  formula  (1),  the  number  of  sensors  that  detect  the  object 
will  be  a  Poisson  random  variable  with  mean 


A 


(3) 


i  r3(-) 

where  the  substitution  A  =  ~  [g(/;i)|'i  obtained  from  equation  (2)  was  made. 


Before  concluding  this  section,  it  should  be  pointed  out  that  the  proposition’s  results 
relating  an  arbitrary  sensor  to  its  nearest-neighbor  sensor  also  apply  to  any  arbitrary  point. 
That  is,  the  distribution  function  Gp(r)  of  the  distance  between  an  arbitrary  point  and  the 
nearest  sensor,  and  the  mean  distance  E(V)  between  an  arbitrary  point  and  the  nearest  sensor 
are,  respectively, 


3.2  PROBLEM  STATEMENT  FOR  R2 


Theorem  2:  Consider  discs  in  two-dimensional  space  with  centers  distributed  according  to  a 
Poisson  distribution  with  mean  A||A||,  where  ||A||  represents  the  area  of  the  set  A.  Suppose 
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that  the  radii  of  all  discs  are  independent  of  the  location  of  the  center  of  the  disc  and  distributed 
according  to  F(r )  with  density  /(r)  and  finite  second  moment.  Then,  the  number  of  discs 
that  cover  a  point  x  is  a  Poisson  random  variable  with  parameter 


Xtt 


(4) 


The  proof  of  this  result  follows  the  proof  of  Theorem  I  provided  in  appendix  B. 

Proposition:  For  the  two-dimensional  case,  the  solutions  to  the  nearest-neighbor  problem  are 

Fd(x)  =  1  —  exp  {— A7rx2}  and  E(D)  =  l/(2\/A) , 

respectively. 
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4.  EXAMPLES 


Example  1:  Suppose  sensors  have  been  distributed  throughout  some  region  in  space  with 
average  distance  between  a  sensor  and  its  nearest  neighbor  equal  to  E(D)  kilometers,  (a) 
What  is  the  probability  that  a  target  in  the  region  is  detected  by  at  least  N  sensors?  ( b )  What 
is  the  probability  that  a  target  in  the  region  is  detected  by  at  least  two  sensors? 

Solution  (a): 


P{X>N)  =  £ 


An  exp  {—An} 


n=N 


n\ 


N- 1 


=  >-£ 


A"  exp{— An} 


71=0 


n! 


where  A  =  /0°°  r3f{r)dr. 


Solution  (6):  Let  X  denote  the  number  of  sensors  that  detect  the  target: 

P{X  >  2}  =  1  -  (P{X  =  0}  +  P{X  =  1}) 

=  1  —  e-A  —  Ae-A 
=  1  —  (1  +  A)e-A 


=  1  - 


(1  + jrTwiiT  r3/(r)dr) 

exP  JT  r3/(r)dr) 


(5) 


Example  2:  (a)  What  average  distance  is  required  between  each  sensor  and  its  nearest  neigh¬ 
bor  in  order  to  have  an  a-%  chance  that  at  least  two  sensors  cover  an  arbitrary  point  in  the 
volume?  ( b )  How  intensely*  should  the  sensors  be  distributed  over  the  volume  in  order  to 
achieve  this  coverage? 

Solution  (a):  Here,  we  want 

P{X  >  2}  =  1  -  (P{X  =  0}  +  P{X  =  1}) 

=  1  -  e~A  -  Ae"A 
=  a. 

*This  refers  to  the  parameter  A:  larger  values  of  A  imply  a  denser  sensor  field. 
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Thus,  one  first  solves  A  -  log(l  +  A)  +  log(l  -  a)  =  0  for  A  (numerically)  and  then  uses 
equation  (3)  to  obtain 


E(D)  = 


(6) 


Having  found  the  answer  to  (a),  one  can  easily  answer  (6):  Solving  (2)  for  A  gives 


A 


r(j) 

367T  [E(D ) 
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5.  COVERAGE  COMPARISON  STUDY 


This  section  provides  a  comparison  of  the  coverage  provided  by  two  systems  of  ran¬ 
domly  distributed  sensors.  Throughout  this  section,  it  is  assumed  that  the  sensors  are  dis¬ 
tributed  over  a  planar  region,  i.e.,  S  C  R2.  Case  1  considers  a  system  in  which  all  sensors 
have  identical  range  sensitivity:  ra  =  3  units.  Thus,  the  range  sensitivity  density  for  this 
system  is  the  dirac-delta  function  5(r  —  i'a).  For  case  2,  the  gamma  density  is  used  to  model 
the  range  sensitivity  over  all  sensors,  that  is, 


f(r) 


1 

I»/3« 


ra-ie-r/0 > 


r  >  0. 


In  the  example  given,  a  =  2.5  and  (3  =  1.2;  thus,  the  mean  range  sensitivity  is  also  a/3  —  3.0. 
The  gamma  density  is  used  to  model  range  sensitivity  to  convey  the  idea  that  some  sensors 
will  be  better  than  others.  Upon  deployment,  some  sensors  will  be  defective  or  give  very  poor 
performance,  whereas  others  will  perform  very  well.  And,  of  course,  several  will  provide 
average  performance  (figure  1 ). 


Figure  1.  Gamma  Density  Function  with  Parameters  a  =  2.5  and  3  =  1.2 
(  The  line  x  =  3  represents  the  idealized  density  with  support  concentrated  at  3.) 
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In  general,  for  /  equal  to  the  gamma  density,  one  has 


A 


C 

r  (a)(3a 

c 

T(a)Pa 


r«+7-l e~r /0 fa 


T{a  +  7)/?a+7- 


From  section  3.2,  it  follows  that  the  mean  parameter  for  the  homogeneous  case  is 

A 7T  J0°°  r2S(r  —  r0)dr  =  \itr2. 


(7) 

A  = 


In  this  simulation,  100  sensors  are  uniformly  distributed  over  a  100  by  100  square 
unit  area.  For  the  homogeneous  system,  with  all  sensors  having  range  up  to  r0  =  3  units, 
the  predicted  coveraget  is  p  =  1  -  exp(— Xth'q)  =  0.246.  On  the  other  hand,  for  the  het¬ 
erogeneous  system  with  range  sensitivity  following  a  gamma  distribution  with  parameters 
ct  =  2.5,/?  =  1.2,  the  predicted  coverage  is  p  =  1.0  —  exp  |  —  r(^„  T(q  -f  2)/?Q+2  j  =  0.326. 
Note  that  the  mean  sensitivity  range  for  the  heterogeneous  sensors  is  a/3  =  3.  Figure  2 
compares  the  coverage  provided  by  just  one  simulation  of  a  randomly  distributed  system  of 
homogeneous  and  heterogeneous  sensors.  The  plot  on  the  left  depicts  100  randomly  placed 
sensors.  The  actual  coverage  provided  by  this  system  is  0.256.  The  plot  on  the  right  contains 
100  sensors  with  0.339  coverage.  This  is,  of  course,  just  one  instantiation  of  the  two  systems. 
In  the  next  section,  data  from  several  instantiations  are  used  to  test  the  long-term  accuracy  of 
the  simulations. 


*Here,  the  parameters  C  and  7  depend  on  the  dimension:  if  the  search  region  is  in  R 2,  then  7  =  2  and  C  =  A7r;  if  the  search  region  is  . 
a  subset  of  R 3,  then  7  =  3  and  C  =  |A7t. 

*The  term  predicted  coverage  denotes  the  probability  that  an  arbitrary  point  in  the  space  is  within  range  of  at  lease  one  sensor. 
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100  Homogeneous  Sensors 


Homogeneous  Sensor  Field 


100  Heterogeneous  Sensors 
Heterogeneous  Sensor  Field 


Actual  coverage  =  .256  Actual  coverage  =  .339 

Figure  2.  Comparison  of  Homogeneous  and  Heterogeneous  Sensor  Fields 


Figure  3  compares  the  number  of  sensors  deployed  to  the  predicted  coverage.  From  the 
formula  for  p,  it  is  clear  that  the  critical  factor  in  the  coverage  statistic  (for  the  case  of  a  planar 
field  of  sensors*)  is  the  mean  coverage  area  of  the  sensors:  r2f(r)dr. 


No  of  Sensors  Deployed  (100  by  100  sq  unit  area) 


Figure  3.  Number  of  Sensors  Deployed  Versus  the  Probability  of  Detection 

by  at  Least  One  Sensor 

'For  a  network  of  sensors  over  a  volume,  the  mean  volume  f£°  r3f(r)dr  is  the  critical  statistic. 
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6.  MONTE  CARLO  STUDY  OF  THE  COVERAGE  STATISTIC 


Suppose  that,  for  each  simulation,  one  randomly  distributes  sensors  over  a  search  region 
and  then  randomly  chooses  a  point  in  this  region.  The  randomly  chosen  point  should  be 
covered  by  at  least  one  sensor  approximately  p  =  1  —  exp  (—A)  percent  of  the  time.*  For  the 
i//i  simulation,  let  =  1  if  the  randomly  chosen  point  is  within  range  of  at  least  one  sensor; 
otherwise,  set  Xt  —  0.  The  random  variables  Xt  are  binomial  with  success  probability  p.  The 
sample  average  P  =  jj  ^2iL\  of  the  simulations  can  be  used  to  determine  the  number  of 
simulations  required  so  that  p  is  within  5%  of  p  with  90%  confidence.  Specifically,  one  wants 
to  estimate  N  such  that 


P(|  p-p  |<  0.05/0)  =  0.90. 


Note  that 


Pp  =  P 

oj  =  p{l~p)/N. 


By  the  Central  Limit  Theorem,  for  N  large,  p  is  approximately  normally  distributed 
with  mean  p  and  variance  p(  1  —  p)/N.  Hence, 


0.90  =  P{\  p-p  |<  0.05p) 

=  p(  -  P~P 
\  v/p(1  ~P)/N 

=  P\\Z\<  0.05p 


<  0.05p 


Vn 


VpU~p), 


y/N 


VpO-  -p) 


where  Z  is  the  standard  normal  random  variable.  In  order  for  p  to  be  within  5%  of  p  (90%  of 
the  time),  it  is  required  that 


0.05p 


t/n_ 

VpO ~-~p) 


1.645; 


that  is,  N  «  (1. 645/0. 05)2(1  -  p)/p. 

^Recall  from  Theorem  1  that  the  number  of  sensors  that  detect  a  randomly  chosen  point  in  space  is  Poisson  with  mean  A.  Hence,  the 
probability  that  the  point  is  not  detected  is  exp  (-A).  The  probability  that  at  least  one  sensor  covers  the  point  is  1  -  exp  (- A). 
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The  analysis  was  performed  for  A  =  0.05  over  a  100  by  100  square  unit  area.  Thus, 
each  simulation  has  roughly  10000A  =  500  sensors.  As  in  the  previous  section,  each  sensor 
in  the  homogeneous  field  has  a  range  sensitivity  of  r0  =  3  units,  and  the  sensor  range  density 
for  the  heterogeneous  system  is  gamma  with  parameters  cv  =  2.5,/?  =  1.2.  The  predicted 
coverages  provided  by  the  two  systems  are  0.757  and  0.862,  respectively.  Note  that,  for  this 
example,  detection  in  the  homogeneous  case  is  a  more  rare  event,  hence  more  simulations  are 
required  to  approximate  the  true  value:  N  =  (1. 645/0. 05)2(1  -  0.757)/0.757  =  348  versus 
N  =  (1. 645/0. 05)2(1  —  0.862)/0.862  =  174  for  the  heterogeneous  case. 


Table  1.  Monte  Carlo  Simulation  of  Coverage  Probability 


System 

No.  Simulations 

P 

P 

Percentage  Difference 

Homogeneous 

348 

0.757 

0.744 

1.72 

Heterogeneous 

174 

0.862 

0.845 

1.97 
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7.  SUMMARY 


Through  an  entirely  probabilistic  approach,  coverage  statistics  for  a  system  of  randomly 
distributed  sensors  having  heterogeneous  range  sensitivities  have  been  derived  in  this  report. 
It  has  been  demonstrated  that  the  case  in  which  all  sensors  have  identical  range  sensitivity 
is  just  one  idealized  case  of  the  heterogeneous  problem.  Simulation  studies  that  support  the 
theoretical  results  are  provided. 
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APPENDIX  A 

CONDITIONAL  UNIFORMITY  AND  INDEPENDENCE 


Lemma:  Let  N  be  a  spatial  Poisson  process  with  parameter  A.  Then,  under  condition 
N(A)  =  k  for  ||A||  >  0,  the  k  points  are  independent  and  uniformly  distributed. 

Proof:  Let  A\,  A2, . . . ,  An  be  n  disjoint  regions  with  U"  A,  =  A.  Suppose  also  that  the 
integers  kx ,  k2,  ■ . .  kn  satisfy  kx  +  k2  +  ...  +  kn  =  k.  Then, 


P{N(AX) 


In  particular. 


kuN(A2)  =k2,...,  N(An)  =  kn  |  N(A)  =  k} 
P\MA x)  =  kj ,  N(A2)  =  fca,  ■ .  ■ .  N(An)  =  kn} 
P{N(A)  =  k} 


[AlMilll 

fc.! 


!1.e-M\Ai\\]  [lAM2|!l!ie-A||/t2|| 


M^nlll*"  -All  An\ 
kn  • 


WMte-\\\A\\ 
fe!  c 


k\ 

[Mill] 

k  i 

[M2III 

k2 

riiA,in 

kx\k2\...kn\ 

Lmii  J 

LmiiJ 

Mil  J 

(A- 1 ) 


P{N{AX)  =  1  |  N{A)  =  1} 


P{N{AX)  =  1, N(A*{)  =  0  |  N(A) 
P{N(Al)  =  l,N(Ac1)  =  0} 


1! 

TTo[ 


P{N(A)  =  1} 


[Mill] 

1 

rMfin 

Mil  J 

L  Mil  J 

1} 


(A-2) 


According  to  (A-2),  given  that  an  event  has  occurred  in  volume  A,  it  is  equally  likely  to  be 
found  anywhere  in  A.  Also,  note  that 


k\ 

(A-3) 

Thus,  equation  (A-l)  implies  independence.  Note  that  the  leading  coefficient  in  (A-l)  is  the 
sum  over  all  the  configurations  in  which  the  k  points  can  be  divided  into  n  subgroups  with  kx 
points  in  the  first  group,  k2  in  the  second,  etc.  Each  of  these  configurations  has  probability 

l|A1f,||A2f2...||A„f"/Pllfc- 


P{N(Ax)  =  kx\N(A)  =  kx}  = 


Mi 
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APPENDIX  B 

PROOF  OF  MAIN  THEOREM 


Theorem  1:  Consider  a  set  of  omnidirectional  sensors  in  three-dimensional  space  scattered 
throughout  some  region  S  C  R3.  Suppose  that  the  sensors  are  spatially  distributed  according 
to  a  Poisson  point  process  with  parameter  A.  Suppose  also  that  the  range  sensitivity  of  all 
sensors  is  distributed  according  to  the  cumulative  distribution  function  F(r)  with  density 
F'(r )  —  f(r)  having  finite  third  moment.  Finally,  assume  that  range  sensitivity  is  independent 
of  the  sensor  location.  Then,  the  number  of  sensors  that  will  detect  a  target  located  at  some 
point  x  e  <S  is  a  Poisson  random  variable  with  mean 


A 


Independence  Assumptions:  ( 1 )  Sensors  are  distributed  according  to  a  Poisson  random  vari¬ 
able.  This  means  that  the  number  of  sensors  in  non-intersecting  regions  of  space  (e.g.,  non¬ 
intersecting  spherical  shells  centered  about  any  point)  are  independent.  It  is  useful  to  think 
of  each  sensor  and  its  sensitivity  extent  as  a  sphere,  with  the  center  of  the  sphere  at  the  sensor 
location,  and  radius  equal  to  the  sensitivity  range  of  the  sensor.  (2)  According  to  the  theorem, 
the  radius  of  each  sphere  is  assumed  to  be  a  random  variable  independent  of  the  location  of 


the  sphere. 


Proof:  Fix  the  origin*  at  any  point  in  R3.  Let  S(r)  be  the  sphere  of  radius  r  with  center 
at  the  origin.  Also,  let  S(r,  r  +  A r)  denote  the  volume  of  the  shell,  or  region,  between  two 
concentric  spheres  centered  at  the  origin  and  having  radii  r  and  r  +  A r,  respectively.  Then, 


S(r,r  +  Ar)  =  ||5(r  +  Ar)||  -  ||S(r)|| 

3 


-7r  [r3  +  3r2Ar  +  3r(Ar)2  +  (Ar)3  —  r3] 

O 

-7rA r  [3r2  +  3rAr  +  (Ar)2] 

O 

47rr2Ar  4-  o(Ar). 


Therefore,  the  probability  of  a  sphere  occurring  in  the  shell  (i.e.,  having  center  in  the  shell 
S(r,r  +  Ar))  with  radius  extending  out  to  the  origin  is  the  product  of  ( i )  the  probability  a 

*For  argument's  sake,  the  origin  is  used,  although  any  fixed  point  would  suffice.  It  might  be  useful  for  the  reader  to  consider  this  point  a 
target,  as  it  will  be  shown  that  the  number  of  sensors  that  detect  this  point  follows  a  Poisson  distribution. 
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spherical  center  occurs  within  the  shell  (i.e.,  A  (47rr2Ar  +  o(Ar)))  and  (ii)  the  probability 
that  the  sphere’s  radius  extends  at  least  out  to  r  (i.e.,  /r°°  f(p)dp): 


That  is, 


P  {N(S(r,  r  +  Ar))  =  1}  =  A(r)Ar  +  o(Ar), 


(B-l ) 


where  N(S(r,r  +  A/  ))  is  the  number  of  sensors  (i.e.,  sphere  centers)  to  occur  in  the  shell 
S(r,  r  +  Ar)  with  radius  extending  out  to  the  origin,  and  A(r)  =  47rAr2  Jj00  f(p)dp. 

It  remains  to  be  shown  that  the  probability  of  the  occurrence  of  two  or  more  spheres 
with  centers  in  the  shell  S(r,  r  +  Ar)  and  radii  extending  out  to  the  origin  is  o(  Ar): 


[47rAr2  Ar  f(p)dp] 2 


<  exp  (— 47rAr2Ar) 


2! 


(B-2) 

(B-3) 


=  o(Ar). 


Note  that  for  small  Ar,  the  exponential  factors  in  (B-2)  are  approximately  1.  The  coefficient 
of  ( Ar)2  in  the  third  factor  is  a  constant. 

From  the  two  independence  assumptions,  it  follows  that  S(r)  has  independent  incre¬ 
ments.  Now,  (B-l),  (B-3)  and  the  independent  increment  property  imply  that  the  number  of 
spheres  with  centers  approximately  a  distance  r  from  the  origin  and  having  radius  out  to  the 
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origin  is  a  Poisson  random  variable  with  mean  A(r).  That  is,  the  number  of  sensors  with  range 
out  to  the  origin  is  a  non-homogeneous  Poisson  random  variable  with  local  intensity  function 
A (r).  Finally,  one  integrates  A (r)  over  all  ranges  to  get  the  cumulative  intensity  function: 


roc  roc 

/  47rAr2  /  f(p)d.pdr 

roo  n 

/  M 

Jo  Jo 

r°°  4 

Jo  3 


4tt\v  dr  dp 


7tV3/(p)^P- 


Since  the  point  chosen  as  origin  was  arbitrary,  the  result  follows  for  any  point  in  Hs. 
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